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Abstract 

The unknown constants in Chiral Perturbation Theory needed for an 
all orders analysis of the polarizabilities and electromagnetic corrections to 
the masses of the pseudo-Goldstone bosons are estimated at leading order 
in l/Nc- We organize the calculation in an l/A'c-expansion and separate 
long- and short-distance physics contributions by introducing an Euclidean 
cut-off. The long-distance part is evaluated using the ENJL model and the 
short-distance part using perturbative QCD and factorization. We obtain 
very good matching between both. 

We then include these estimates in a full Chiral Perturbation Theory 
calculation to order for the masses and for the polarizabilities. 

For the electromagnetic corrections to the masses, we confirm a large vio- 
lation of Dashen's theorem getting a more precise value for this violation. 
We make comparison with earlier related work. Some phenomenological 
consequences are discussed too. 
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1 Introduction 



Virtual electromagnetic (EM) effects in purely strong processes can be important 
in precision situations. This is especially true in the case of isospin breaking 
contributions to hadron masses and some hadronic processes. If we want a high 
precision description of the latter, we need to know not only the effects due to the 
quark mass difference — m„, which is a quantity we would also like to extract 
from these experiments, but the size of the electromagnetic contributions. 

That these contributions can be sizeable in certain cases is best illustrated 
in the case of the observed tt^ — 7r° mass difference which is almost entirely due 
to photon loops [|I|]. At present, we cannot directly use QCD to estimate these 
effects. Some first progress using lattice QCD has been made recently in 0. 
Instead we turn to the method of Chiral Perturbation Theory (CHPT). In the 
purely mesonic sector for the strong and semi-leptonic processes, this was started 
as a systematic program by Gasser and Leutwyler || and has since then been 
extended to a large variety of processesp. In the case of 77 — ^ ttttQ, 0|, nn- 
scattering [§] and vr iu^y this has even been done to the two-loop order. 
In TTTT-scattering, electromagnetic corrections might also become relevant at the 
precision achieved by the two- loop order calculation. They are already important, 
at the present level of precision, in various others form factors. 

Chiral Perturbation Theory for virtual electromagnetic effects was first de- 
scribed at lowest order (e^p°) in |T^. Urech[^] has recently systematically stud- 
ied the next-to-leading order terms. His work has been mainly dedicated to the 
EM correction in the masses of the lowest pseudoscalar Goldstone bosons. This 
program has been later expanded into a few more form factors by Neufeld and 
Rupertsberger [T3|| . The results were, however, quite dependent on the values 
for the new coupling constants appearing at order in the chiral expansion. 
Unfortunately, contrary to the case of strong and semileptonic processes, it is 
impossible to determine all these constants from experimental data. The calcula- 
tion of these constants is the main subject of this paper. In addition we provide 
estimates for the counterterms appearing in CHPT at order in the processes 
77 — > PP. Here we extend the work done in [|1^] to the charged case and make 
predictions for the polarizabilities of pions and kaons too. This side aspect is 
discussed in Section |^. 

The real problem in calculating the constants is that it requires an integra- 
tion over internal photon momenta. For instance, in the estimate done for the 
strong sector in |jl5[ using the lowest lying resonance saturation, only their lowest 
order couplings are relevant. When one integrates over all photon momenta, the 
couplings of the resonances to all orders need to be known, thus making these 
estimates much more difficult. The lowest order constant, called C in Section 
1^, has a very long history, it was first estimated in [|^ using PCAC and then 
saturating two-point functions with resonance exchanges. A different technique 
was subsequently developed by Bardeen, Buras and Gerard in ||16i for weak non- 
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leptonic matrix elements, now generally referred to as the l/Ai'c-approach. This 
approach was then used together with the Das et al. sum rule to estimate the 
pion mass difference at lowest order, or equivalently C, in |T^. Here a proper 
separation of long and short distance contributions was also possible. The same 
method has been used in the chiral quark model |18[ and in the extended Nambu- 
Jona-Lasinio (ENJL) modelfl^. All of these only allowed estimates of the lowest 
order constant C since they were based on the Das et al. sum rule. 

The calculation of the pion and kaon EM mass differences beyond lowest order 
in CHPT, using saturation by the lowest lying resonances, was recently performed 
by Donoghue, Holstein, and Wyler in pO[ and Baur and Urech in |21|. For earlier 



attempts see . In [g^ the chiral logs aXv=l GeV were used to estimate these 
EM mass differences. These papers all had to make assumptions about the short 
distance behaviour. More comments about these assumptions are in Section ^ 
The short-distance contribution was introduced in an operator-product expansion 
(OPE) framework in p4|. This is discussed in Section |0| . 



The method used in this paper is an extension of the original l/N^ method[|T^. 
We use an off-shell Green function. This method was used by us previously in 
the calculation of the AS* = 2 hadronic matrix element in the K^-K^ system and 
commonly parametrized by the so-called factor ||25|| . 

The paper is organized as follows: in Section ^ we shortly discuss CHPT for 
electromagnetic corrections and define our set of counterterms. Our set is some- 
what more appropriate to the large Nc limit and different number of flavours. We 
also point out in some detail the gauge dependence of the generating functional. 
No observable quantities do of course depend on the gauge but the infinite parts 
of the constants at next-to-leading order do depend on it. In Section ^ we explain 
the method. This section also gives a short overview of the ENJL model that 
we use for the long-distance contributions. The main contributions are the latter 
due to the photon propagator. In Section || we give the results and use CHPT 
at leading order in l/N^ to extract the CHPT constants. We compare with the 
earlier work in Section ^. Section |^ summarizes the consequences for the ratios of 
the light current quark masses. Section |^ contains the results on pion and kaon 
polarizabilities and we present our main conclusions in Section M. 



2 Chiral Perturbation Theory Analysis 

In this Section we use CHPT to analyse the two-point functions, 

n(g2) = zjd^x e'"" (0|T [P{0)P\x)) |0) , (2.1) 

in the presence of electromagnetic interactions. We shall do this to first order in 
aQED and to order p"^ in CHPT. The pseudoscalar source P{x) in ( |2.1D is defined 
as ijiJ'lbQb) with a, h quark flavour indices and colour indices summed inside 
the parenthesis. 
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At lowest order in the chiral expansion 0{p'^) the strong interactions 
between the lowest pseudoscalar mesons including external vector, axial-vector, 
scalar and pseudoscalar sources are described by the following effective Lagrangian 



C 
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{tr (d^UD^'U^) + tr (xU^ + Ux^)] 
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where denotes the covariant derivative 
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and a^{x) 



are external 3x3 vector and axial- vector field ma- 

0. 



In Eq. 

trices. When electromagnetism is switched on, v^{x) = \e\QAfj_{x) and 
Here A^{x) is the photon field and the light-quark electric charges in units of the 
electron charge |e| are collected in the 3x3 flavour matrix Q = |diag(2, —1, —1). 
In Eq. ( p.2|) x = 2_Bo(A1 + s{x) + ip{x)) with s{x) and p{x) external scalar 
and pseudoscalar 3x3 field matrices and Ai the 3x3 flavour matrix Ai = 
diag(mtj, m^, m^) collecting the light-quark current masses. The constant Bo is 
related to the vacuum expectation value 



\m\o)l 



=u,d,s 



-F^Boil + Oim,)) . 



(2.5) 



In this normalization, Fq is the chiral limit value corresponding to the pion 
decay coupling ~ 92.4 MeV. In the absence of the U(1)a anomaly (large 
limit) [^, the SU(3) singlet rji field becomes the ninth Goldstone boson which 
is incorporated in the $(x) field as 



To this order, the two-point functions in (|2.1| ) have the following form 

9 r2 p2 

n(g^) = (2.7) 



where niQ is the pseudo-Goldstone boson masses to that order, i.e. ml = Bo{ma + 
rUb), with rria the flavour a quark mass. Here we are interested in the isospin 
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breaking corrections to the poles of these two-point functions induced by electro- 
magnetism to order e^. So we will set rriu = in the calculation. 

To take into account virtual photons excitations we need to add to the La- 
grangian in ( p.2|) the corresponding kinetic and gauge fixing terms, i.e. 

- If'^'^F,^ - ^(Y^ (d.A'f (2-8) 

where F'^'^ = d^A^ — d'^ is the electromagnetic field strength tensor. The 
parameter ^ is the gauge fixing parameter which is zero in the Feynman gauge, 
one in the Landau gauge and four in Yennie's one. 



2.1 Lowest Order Contribution 

At lowest order in the chiral expansion, electromagnetic (EM) virtual interactions 
of order between the pseudo-Goldstone bosons are described by the following 



effective Lagrangian P, 27, 28 



= e^C^{^l)ti(Q^)+e^C2{^l)ii(QUQU^). (2.9) 

Here the U field matrix is the U(3) symmetric one in (|2.6| ) and $o = & + Fq, 
where 9 is the so-called QCD theta- vacuum parameter. The constant C = 6*2(0) is 
the coupling introduced in ||15| . There are no loop contributions to this order and 
therefore Ci^2(0) and the derivatives of Ci^2{x) at x = are finite counterterms 
not fixed by symmetry alone. In the large Nc limit, Ci,2(0) are of order Nc whereas 
the n-th derivative of Ci^2{x) is of order 1/N^"'~^. To order e'^p^, the correction 
to the pole position of the two-point function in ( ^.7| ) is zero for P = tt^ and 
P = rjQ, while charged pion and kaons get the same non-zero correction, namely 



2C 

mlM{n+,K+) = T^s • (2-10) 



This is the so-called Dashen's theorem [rU|. The pole position for the rji-rji two- 
point function gets corrected by 

_ sjcm+cm) .2 11) 

from EM virtual interactions. 



2.2 Next-to-Leading Order Contribution 

In this Section we shall report on the CHPT order p'^ virtual EM corrections. 
The first type of these corrections is the emission and absorption of a photon 
by a pseudo-Goldstone boson line (see Fig. |T]). To the order we are interested 
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Figure 1: The photon loop contributions to (|3.1|) . The crosses are pseudoscalar 
currents. The wiggly line is the photon. The full lines are pseudoscalars. 



here, the '~)P^P~ and 'yyP^P^ vertices come from the Lagrangian in ( |2.2| ). This 
contribution needs a counterterm of order e^p^ to make it UV finite. 

The complete CHPT order effective Lagrangian describing virtual EM 
interactions of order between pseudo-Goldstone bosons and external vector, 
axial- vector, scalar and pseudoscalar sources was written in |TT[ (see also |P^). 
The coefficients of this Lagrangian are the needed counterterms to absorb the UV 
divergences of order e^p^. To construct this Lagrangian, some Cayley-Hamilton 
relations for SU(3) matrices were used. We would like here to include also the 
ninth pseudo-Goldstone boson as above and work with U(3) matrices (see Eq. 
( |2.6| )), which is the symmetry in the large Nc limit (A^^^ is the number of colours). 
This is useful for our calculation since we want to use the 1/iVc-expansion as the 
organizing scheme. To order e^p^, one has to add to the Lagrangians in 
( ^TSf ), and (^I9D the following one: 









[d^'u^d^u) tr (g2) 
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e^K2 




[d^'UW^u) tr (QUQU^) 
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e'Ks 
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:^l)F^tr 
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:^l)F^tr 


[d^'U^D^UQU^QU + D^'UD^U^QUQU^) 


+ 


e'K, 


:^l)F^^T 


[x^U + U^x) tr (Q') 
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e^Ks 




[x^U + U^x) tr (QUQU^) 


+ 


e^Kg 




[{x^u + u^x + xu^ + ux^) Q') 


+ 




)(*o)^otr 


{{x^u + u^x) QU^QU + [xu^ + Ux^) ququ^) 


+ 


e^Kn 


{^l)F^ ti {{x^U - U^x) QU^QU + {xU^ - Ux^) QUQU^) 
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+ e^K^2i^l)F^tr(UD^U^[v^Q,Q] +U^D^'U 

+ e'K,,i^l)F^ tr (vjQV^^g + V^QV^Q 

+ e^K,^{^l)F^tT ip^UQU^Q - D^QUQ 

f tr {QUQU^ 



+ e'Kn{<^l)F^(ti{Q 
+ e'k,smFf(tr{Q 



+ e'K,,i^l)F^{tT{Q' 

+ e2$oi^2o(*o)^otr 

+ e2$oK23($o)^otr 

+ e2$oK24($o)^otr 



+ tr [QUQU^ 
+ tr {QUQU^ 
f/M tr fg2 



X^U-U^x) tr {QUQU^ 



(x^U + f/tx) QU^QU - (xU^ + f/tx) QUQU^ 
X^U-U^x) QU^QU~ (xU^ 



Ux^) QUQU\ 

{2.12) 

It is interesting to study the large behaviour of the coefficients in ( |2.12|) . In 
the large A^'^ limit, the couplings i^i(O), i^2(0), ^^7(0), K^ifi), i^i5(0), Ki6(0), 
-^'22(0), -^'23(0), and ^^24(0) are order 1/Nc and the couplings -^'17(0), Ki^{<S), 
-^'19(0), -^'20(0), and -^'21(0) are order l/N^- The rest of the constants Ki{0) and 
Ki{0) in ( |2.12| ) are of order 1. Each derivative of the Ki{x) and Ki{x) functions 
brings in an additional factor In the rest of the paper, we call the functions 

Ki{^l) and Kj{^l) at $0 = the couplings Ki and Kj. 
The covariant derivatives are defined as follows, 



and the symbol 6^0 is 



6 J = dj + 2tr (a. 



(2.13) 



(2.14) 



The relation between these couplings and the ones defined in [|rT| |l^ is 



A', = A'l + A's 

K, = -k; ; K, = K, 
K, = k,-2ks ; Ke = ke~k^. 

From here one gets the following large Nc limit relations 

Ki = -K3 and 2K2 = K4. 



(2.15) 



(2.16) 
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The rest of the tilded couphngs introduced here were not included in those ref- 
erences because they worked in an octet symmetry framework. The expressions 
of the corrections to the poles of pseudoscalar two-point functions at this order 
are given in Appendix 0. The coefficients of the effective Lagrangian in (|2.12| ) 



are again not determined from chiral symmetry arguments alone. Its estimation 
is the central subject of this work. We explain the technique we use in the next 
section. 

2.3 Gauge Dependence of the Various Quantities 

A subtle issue is involved here. The generating functional in terms of colourless 
external fields, as used in , is not independent of the gauge chosen for the gauge 
fields propagators (in particular the photon one). As a well known consequence 
Green functions are not gauge invariant in general. The underlying reason is 
simply that external sources are in general charged so they transform under the 
gauge group non-trivially. Of course, for observable quantities, like the mass shift 
we obtain from the two-point function or any other physical quantity, the 

result has to be gauge invariant. So this gauge dependence disappears when the 
external sources are on the mass-shell. 

For instance, one obtains a gauge dependence in the result for the two-point 
function which only cancels when the meson created by P^(x) and destroyed 



by P{x) is on the mass-shell. This means that some of the couplings Ki and Ki 
are actually U(l) gauge dependent. In practice, in our CHPT calculation we 
fix the gauge for the photon propagator to be the Feynman one (i.e. ^ = 0). 
The same gauge was used in the CHPT calculations in |Tl|, |T^. If one wants to 
compare or use the values of the couplings we get with the ones obtained from 
experiment or other model calculations, one should make sure that the Feynman 
gauge (the gauge we used) is used in the CHPT calculation or in the model 
calculation. Alternatively, one can of course compare directly the same physical 
quantity. 

This gauge dependence does not reduce the number of parameters in the 
CHPT Lagrangian, since choosing a clever gauge fixing in order to remove a 
constant from the Lagrangian would bring back the parameter in the photon 
propagator. 



3 Calculation of the Counterterms 



We calculate directly the two-point function defined in ( p.ip in the presence of 
EM interactions to order e^. In practice, this means the calculation of 
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Where Ilp'pyy{q,r) is the following four-point function 

U^ji'pyy{q,r)=t^ J d^x J d% J d^2e-'(''"+"(^-")) {0\T [P{0)P\x)V^{y)V''{z)) |0) 

(3.2) 

Here V^(x) = {q{x)Q'y^q{x)) with q{x) the SU(3) flavour vector {u{x), d{x), s{x)). 
Notice that since the photon momentum r is integrated out, we should know this 
function at all energies. In particular, this means to all orders in a low energy 
CHPT expansion of this function in the momentum r. In order to extract as 
much information as possible of the two-point function ( |3.1| ), we calculate it at 
off-shell values of as well. 

Let us now discuss on the U(l) gauge invariance of ( |3.2| ). When the external 



pseudo-Goldstone bosons are on-shell, the four-point function in (|3.2| ) is U(l) 
gauge covariant and fulfills 

rp.'^p'pvvil^ = ruTi-p'pvviQ^ = 0. (3.3) 



Therefore the gauge dependent term proportional to ^ in ( p.l|) cancels. This is no 
longer true when we move to off-shell values. In that case, the ^ term gives a 
non-zero contribution since we have no gauge covariance for ( |3.2| ) (see comments 
in Section pl3|) . 



The technique we use here is similar to the one introduced in W& and used in 



[|r^ , and is the variant presented in |2^ . This consists in introducing an Euclidean 



cut-off in the integrated out photon momentum. This cut-off /i both serves to 
separate long and short distance contributions and as a matching variable. After 
reducing the two pseudoscalar legs of 11 (g^) in ( |2.1| ) the result has no anomalous 
dimensions in QCD. We have then to find a plateau in the cut-off /i if there is 
good matching. The photon propagator will help to produce it. 

So, after passing to Euclidean space Eq. (p.l|), we introduce the cut-off fj, in 
the photon momentum r^; as follows 

1 pfj, poo 

drE = drE + drE- (3.4) 
Jo J^i 

For the short- distance part we can perform the full calculation, see next Section. 
In particular, we have obtained, in the large Nc limit and up to order 
the short- distance contributions to all the terms in (|2.12|) and not only those 
accessible via (B.!!). 



3.1 Short-Distance Contribution 



The higher part of the integral in ( |3.4| ) collects the contributions of the higher 
than /i modes of the virtual photons. The effective action obtained by integrating 
out the virtual photons with modes higher than in ( |3.1| ) can be expanded in 
powers of like in an OPE. We compute these contributions up to order 
This contribution was first introduced in . 
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(a) 



(b) 





(c) 



Figure 2: The short distance contributions. The wiggly hne is the photon, the 
curly one the gluon, the full line a quark. There are also the crossed and sym- 
metric configurations. 

There are four types of contributions. There is a pure QED fermion mass 
renormalization contribution (see Figure 0(a)) due to the fact that the quark mass 
in QED runs proportionally to the mass itself. The QED log divergency produces 
a contribution to the effective action proportional to ln(A/ /i) where A is the scale 
where the input current quark masses are renormalized in QED. Within the ENJL 
model the input current quark masses are encoded in the values of the constituent 
quark masses. We have fixed the physical values for the constituent quark masses 
by comparing the ENJL predictions to low energy observables typically at some 
scale between the rho meson mass and the ENJL model cut-off Aenjl = 1-16 
GeV. Accordingly, we will vary the scale A in that range. Of course, there will 
be some kind of double counting since we cannot disentangle the EM virtual 
contributions to the ENJL parameters from the experimental values. But they 
give order e'^ corrections. The contribution of the QED self-energy diagram in 
Figure ^(a) to the effective QCD Lagrangian, in the presence of EM virtual 
interactions, is 



3 a 




(3.5) 



2 vr 
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where colour indices are summed inside parenthesis. Its effective reahzation can 
be calculated to all orders in an l/N^. expansion since it can be written in terms 
of just bilinear QCD currents. At low energies, in terms of the lowest pseudo- 
Goldstone bosons and external sources and at lowest order in the chiral expansion, 
it only contributes to the couphngs Kiq and Ku in ( p.l2| ). 



Kfr = -K?r = ^,in[-]. (3.6) 




Notice that here enters the scale A where the current quark masses are defined. 
The remaining part of the integration from A to cxd is absorbed in the definition 
of the current quark masses. Here we can also indicate the type of corrections 
existing to (|3^). First, at the quark-gluon level there are and aas correc- 
tions to (|3.5|) . Then when going from the quark-gluon level expression in (|3.5| ) to 
the hadronic level one, there are corrections that are higher order in chiral power 
counting. To obtain ( p.6|) we used the order strong Lagrangian in (|2.2|) . There 
are thus CHPT 0{p^) corrections to this result. 

The presence of explicit dependence on A in the short-distance contribution 
to KiQ and Ku indicates that one has to be careful when using rules of the strong 
sector to obtain naive estimates of order of magnitude sizes of CHPT parameters. 
This we will refer to later as failure of naive power counting. 

This problem will appear whenever non-leptonic couplings of interactions 
other than the strong interaction come into play. In particular it also shows 
up in weak non-leptonic decays. There the effects are suppressed by extra inverse 
powers of the ly-boson mass, so its numerical importance is negligible. 

A similar contribution comes from diagram (b) in Fig. ^. Here it is not 
the scalar and pseudoscalar current that are renormalized but the vector and 
axial-vector ones. They have to be defined at the scale A and again the QED 
running can lead to a contribution. The cross in Fig. ^(b) denotes an insertion 
of an external vector or axial- vector current. This will contribute to K12 and, 
with one more insertion, to Ki^ and K14. There are in principle short-distance 
contributions of this type but they vanish because of the global chiral invariance in 
QCD perturbation theory. There will be short-distance contributions from this 
diagram but to higher order operators in CHPT, for instance to magnetic-like 
structures, etc. So we have 

k!? = k!^ = k!^ = o. (3.7) 

There can be long-distance contributions of this type due to the spontaneous 
breaking of the axial symmetry. 

The third and fourth type of contributions to the effective Lagrangian are 
of order N^as- These are the well known box- and penguin-type diagrams (see 
Figs. H(c) and ||(d) and crossed versions). To order they contribute to the 
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effective Lagrangian as follows. The box-type contribution is 

a,b,c,d 

+ 4 (rLQ^VR) (rRQ-'ql)] . (3.8) 
and the penguin-like contribution is 

a,b,c 

- 2 (qlqR) {qRiQTQl) - 2 (r^^^) (^l(Q')'^?S)] • (3.9) 

Latin indices are for the quark flavours and colour indices are summed inside the 
parenthesis. The meaning of ag^{fi^) is given at the end of this section. Penguin 
diagrams with photon and gluon interchanged does also exist but there the photon 
is at low energies and is thus included in the low energy part calculated in Section 



3.2. 



In addition there are also diagrams obtained from inserting external f^, a^, s 
and p currents on the internal quark lines of the box- and penguin-like diagrams 
in Fig. ^. These do not contribute to the effective action at order e^p^ we consider 
here. 

The factorizable contribution (i.e. the leading order in 1/Nc contribution) 
to the realization of these effective Lagrangian, (|378|) and ( p.9|) , can be easily 
obtained since, in that limit, this Lagrangian is just products of QCD quark 
currents. The low-energy realization of QCD quark currents in terms of pseudo- 
Goldstone bosons and external sources can be worked out in terms of the cou- 
plings appearing in the QCD chiral Lagrangian. To order p"^ this chiral La- 
grangian is in ( p.2|) , at order p^, it can be found in [Q. 



Therefore, to order p^ in the chiral expansion and leading order in 1/Nc, the 
contributions to the couplings in (|2.9|) and ( p.l2| ), from integrating out virtual 
photons with Euclidean momentum larger than fi in ( p.l[ ) are 

3<(/.2)FXeff(/^^). 



c 



IT fi"^ 
fi{if,,2\ f2 



3 af{^,') 



9 ^ ^ 2 vr 



6 TT /i 



2 



^ ^QED 9 ,^QED , 3 <(/x^) {2Ls-H2)Bls{f^') 

Kio = ilio -77^9 > All = All +7 2 • 

(3.10) 
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All others are zero because of l/N^ counting. Here, Fq ~ (86 ± 10) MeV is the 
CHPT chiral limit value of the pion decay constant F^^. The constants Lj and 
Hj are couplings of the order strong chiral Lagrangian in [Q. They reabsorb 
the UV divergences that appear in CHPT calculations to order p^. In particular, 
the values of the renormalized finite parts at a CHPT scale z/ = Mp of the 
constants we need in the MS scheme used in 0] are 0: L4 = — (0.3±0.5) x 10"'^, 
L5 = (1.4 ± 0.5) X 10-3, Lg = -(0.2 ± 0.3) x lO'^, Lg = (0.9 ± 0.3) x lO^^ and 
2L8 - H2 = (2.9 ± 1.0) X 10-3 |9[. The scale dependence of the U couphngs 
cancels out when the next order in l/N^. compared to the one in our calculation 
is included. We have chosen the scale z> to be a typical hadronic mass scale. The 
scales z/, v and are unrelated. The final error assigned to the short-distance 
contribution takes into account this 1/Nc uncertainty. 

The coupling does not receive perturbative contributions in QCD and 
QED because it modulates a L'^R^ two-point function. It can however receive 
non-perturbative contributions at leading order in l/Nc- In the case of the ENJL 
model, these appear proportionally to constituent quark masses. 

The expressions in Eq. ( |3.10|) have both order and corrections. The 
product af'^(/i^) i?Qgg(/i^) has to be understood as 

r2\ \ -9/11 



(3.11) 



and 



We use the one-loop large expression 



afifi')^ Tdx^^^^. (3.12) 



2\ l^TT 

= n An f 2/A2 T 3.13 

with Aqcd = 300 MeV and Bo{lGeV^) =1.6 GeV [||]. This corresponds to using 
the large Nc renormalization group to improve the purely perturbative result. 

The short-distance contributions to pseudo-Goldstone boson masses can be 
found in Appendix p. 

We can get an estimate of the contributions suppressed in 1/A^c by keeping 
the leading in 1/Nc contributions from factorization in ( |3.8| ) and ( p.9| ) at the 
quark-gluon level language but including the Lj that vanish for ^ 00 for the 
hadronic realization. This gives an additional set of nonzero terms: 

^ 9 ^ 3 TT ' 



^^ = -9^^ = ' (3.14) 
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3.2 Long-Distance Contribution and Matching 



We need here the two-point function ( |3.1|) for external energies up to the Eu- 
chdean scale fi. This scale is expected to be around 1 GeV, thus beyond the 
applicability of CHPT. We therefore need to resort to models. We have chosen 
the ENJL model for the reasons given below in Sect, 
contribution to the two-point function in 
same ENJL model we use here in 



3.2.1. The low-energy 



was calculated already within the 
I4| for neutral pions. Here we needed to ex- 



tend it to any flavour structure. Since the technique has already been explained 
several times p5|, |30[| we will not repeat it here. 



Afterwards, we integrate over the photon momentum in the Euclidean space 
up to r^; = /i. This gives us the lower part of the integral in 



3.2.1 The ENJL Model 



For recent comprehensive reviews on the NJL |^T| and the ENJL models 



see 



\p3i , |34| . Here, we will only summarize the main features, notation and reasons 
why we have chosen this model. More details and some motivations on the version 
of the ENJL model we are using can be found in |T9|, ^ . 
The kinetic part of the Lagrangian is given by 



C 



J^EN.TL 
kin 



q {i-f'' {df, - iv^ - ia^^jr, - iG^) - {M + s - ijr)^)} q . 



(3.15) 



Here summation over colour degrees of freedom is understood and we have used 
the following short-hand notation: q = (^u, t>^, a^, s and p are external 

vector, axial- vector, scalar and pseudoscalar field matrix sources; M. is the quark- 
mass matrix; is an external colour source transforming as the gluons do in 
QCD. The ENJL model we are using corresponds to the following Lagrangian 



ENJL 



c 



Aenjl 
kin 



a,b 



-9v E [{TlI^qI) {qll^qi) + (g^T^^^y (??j7.?S)] • (3-16) 

a,b 



Here a, h are flavour indices, '^r(l) = (1/2) (1 + (— )75) ^ and 



9v 



STT^Gy (Aenjl) 

JVcJi-ENJL 



9s 



47r^Gs(AENjL) 

JVcJi-ENJL 



(3.17) 



The couplings Gs'(Aenjl) and ^^(Aenjl) are dimensionless and Oil) in the 1/Nc 
expansion and summation over colours between brackets in ( |3.16| ) is understood. 

This model has three parameters plus the current light quark masses. The first 
three parameters are Gs, Gy and the physical cut-off Aenjl of the regularization 
that we chose to be proper-time. Although this regulator breaks in general the 
Ward identities we impose them by adding the necessary counterterms. The 
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light quark masses in Ai are fixed in order to obtain the physical pion and kaon 
masses in the poles of the pseudoscalar two-point functions The values of 
the other parameters are fixed from the results of the fit to low energy effective 
chiral Lagrangians obtained in They are Gs — 1.216, Gy — 1.263, and 



Aenjl = 1.16 GeV from Fit 1 in that reference. Solving the gap equation, we 
then obtain the constituent quark masses: M„ = = 275 MeV and Mg = 427 
MeV. 

The model in (|3.16|) is very economical capturing in a simple fashion a lot 
of the observed low and intermediate energy phenomenology. It has also a few 
theoretical advantages. 

1. The model in ( p.l6| ) has the same symmetry structure as the QCD action 
at leading order in 1/Nc In the chiral limit and for G5 > 1 this 
model breaks chiral symmetry spontaneously via the expectation value of 
the scalar quark-antiquark one-point function (quark condensate). 

2. It has very few free parameters. These are unambiguously determined 
from low energy physics involving only pseudo-Goldstone bosons degrees 



of freedom 35 . 



3. It only contains constituent quarks. Therefore, all the contributions to 
a given process are uniquely distinguished using only constituent quark 
diagrams. Within this model there is thus no possible double counting. In 
particular, the constituent quark-loop contribution and what would be the 
equivalent of the meson loop contributions in this model, are of different 
order in the l/N^ counting. As described in this model includes the 
constituent-quark loop model as a specific limit. 

4. Resummation of fermion-loops automatically produces a pole in all main 
spin-isospin channels within the purely constituent quark picture. This is 
qualitatively the same as in the observed hadronic spectrum. 

5. It provides a reasonable description of vector and axial vector meson phe- 
nomenology [0. 

6. Some of the short distance behaviour is even the same as in QCD. For 



instance, the Weinberg sum rules [3S] are satisfied. These are required in 



some cases for convergency, for instance in the Das et al. sum rule |1|]. 

7. The major drawback of the ENJL model is the lack of a confinement mech- 
anism. Although one can always introduce an ad-hoc confining potential 
doing the job. We smear the consequences of this drawback by working 
with internal and external momenta always Euclidean. 
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3.2.2 Extraction of the Long-Distance Contributions 

The integral over the Euchdean photon momentum is done at fixed r|;, i.e. we 
perform the angular integration first by a Gaussian procedure. We then fit the 
result obtained for the two-point function (|3.1| ) for a fixed value of the external 
momentum to a series of Chebyshev polynomials in r^;. The remaining integral 
over r\. is then straightforward. We have choosen the points in the Euclidean 
region and near = where we expect the artifacts of the model (constituent 
quark on-shell effects) to be suppressed. 

The lower part of the integral in (|3.4|) has a non-analytical component too, 
we have however checked that, in the region of interest for \q^\ < 0.2 GeV^ and 
fE^^ GeV the Chebyshev polynomials give a good fit and the non-analytical 
behaviour is actually very smooth. We have checked this by calculating the lower 
part of the integral within lowest order CHPT with an explicit Euclidean cut-off 
fi. For low fi this agreed with the ENJL calculation. 

The two pseudoscalar legs of the resulting two-point function are then reduced 



off-shell, for details of the reduction technique see |3g]. Since we are working in 
the large Nc limit, this has been done for the fiavour structures uu, dd, ss, ud, 
us and sd. The CHPT expressions for these fiavour combinations in the large N,, 
limit after the reduction are given in Appendix 

As an example of the quality we have shown in Fig. H, the reduced two-point 
function in the chiral limit as a function of for the integral in (|3.4| ) up to = 
0.5 GeV for the charged case. The curvature is purely due to the chiral logarithm. 
Notice that this is well reproduced by the ENJL calculation. Similar good fits 
were obtained for all the other combinations. From the analytic part of the fit 
we can then, for each fiavour case, extract the combination of coupling constants 
as given in Appendix 

3.2.3 Matching 



Summing the short- distance calculated in the Section |3.1| and the long-distance 
obtained in the ENJL model we get the two-point functions in ( |2.lD with the two 
pseudoscalar sources reduced. 

We have studied the matching of the long- and short-distance contributions 
by looking at the stability in the /i scale of the the results. We have plotted the 
charged pion mass difference in Figure ^. The matching is quite good above ener- 
gies around (0.6 ~ 0.8) GeV due to the presence of the photon propagator. This 
is because the presence of the photon propagator is enough to cut the high 
energy contributions. This happens despite the fact that the vector and axial- 
vector propagators obtained within the ENJL have only an acceptable behaviour 
up to around (0.5 ~ 0.6) GeV in the kaon channel. We have checked that this 
stability (matching) region can be enlarged just by imposing the correct high- 
energy behaviour in the vector two-point functions coupling to the photon in this 
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Reduced two-point function 

-0.0022 I 1 1 1 1 r 




-0.0036 I 1 1 1 1 1 1 1 1 

-0.04 -0.035 -0.03 -0.025 -0.02 -0.015 -0.01 -0.005 

q"2 GeV«2 



Figure 3: The long distance contribution to the reduced two-point function as a 
function of for the integral in (|3.4| ) up to /i = r^; = 0.5 GeV. The curve is the 
CHPT fit. All curvature is due to the chiral log. 
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Figure 4: The long-distance (LD), short-distance (SD) and the sum of the 
Nf. ^ oo results as a function of the matching scale /i for m^+|EM- 



casea. The predictions for the couplings remain however mostly unchanged and 
within the quoted errors. This is because the presence of the photon propagator 
suppresses these contributions. 



4 Results 

In this section we give the results of our calculation. In the large A^^^ limit and 
rUu = TTtct there are four independent combinations of pseudoscalar two-point 
functions. We have calculated the following combinations at off-shell values of 
for our analysis: 

1. Neutral case : zero quark masses, equal quark masses corresponding to the 
kaon mass and different quark masses corresponding to the kaon mass. 

2. Charged case : zero quark masses, equal quark masses corresponding to 
the pion mass, equal quark masses corresponding to the kaon mass and 
different quark masses corresponding to the kaon mass. 

^For some work where the correct high-energy behaviour is imposed using QCD-hadron 
duaUty see pO| . 
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Particle 


2 


LD 


SD 


Stability 




10"3 GeV2 


10-3 GeV^ 


10-3 GeV^ 


/i in GeV 


7r+ 


0.95 


0.67 


0.28 


0.60 - 0.85 


K+ 


1.93 


1.47 


0.46 


0.65 - 0.90 




-0.01 


-0.006 


-0.004 


0.65 - 0.80 




0.98 


0.79 


0.19 


0.65 - 1.00 



Table 1: The oo results for the electromagnetic contributions to the meson 

masses squared and the violation of Dashen's theorem. Also quoted are the long- 
and short- distance contribution and the stability region. 

To order e^p^ there appear five combinations of coupling constants. In the long 
distance case we have worked in the ENJL model to all orders in the chiral 
expansion. We will thus also obtain an estimate of the e^p^ and higher corrections 
to the long-distance contributions. From this analysis we have got four of the 
counterterms of the e^p^ Lagrangian in |TT| . Remember there are ten in the large 



Nc limit, where three of them involve external vector or axial- vector sources. This 
is presented in Section TTl . We first extract the relevant corrections to the masses 
directly. 

4.1 EM corrections to the Masses and Dashen's Theo- 
rem 

We take the formulas of Appendix ^ and fit them to the reduced two-point 
function of the relevant particle at a fixed value of /z. We then use the same 
chiral formulas to extrapolate it to the pole. That way we obtain the long distance 
contribution to the various masses. For the short distance we take the results of 



Section p.l| with L4 = Lg = to stay in the Nc ^ 00 limit. The results we obtain 
are given in Table |T] where we also quote the stability region. The contributions 
for the neutral pion are are always very small. The contributions of short and 
long distance are both of course /^-dependent. The numbers given are for the 
middle of the stability region. We have also quoted the result directly for the 
violation of Dashen's theorem given by 

AM|a, = (m^+ - ml, - ml+ + m^o)^^^ . (4.1) 

As an example of the stability we have plotted the long-distance, the short- 
distance and the sum of the contributions to ^^+Iem -^^S- § ^iid similar for 
AM|;jy^ in Fig. ^. The matching for the other quantities is not quite as good but 
quite acceptable. There are several reasons for the very good matching of AM^^j. 
First is that in this combination, the leading effect the QED quark mass renor- 
malization only appears multiplied by pion masses. In addition, though terms 
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Figure 5: The long-distance (LD), short-distance (SD) and the sum of the 
Nf. ^ oo results as a function of the matching scale /i for AM|j^. 

like m^- appear for individual pseudo-Goldstone boson EM corrections (in the 
charged pion mass for instance), they drop out in the combination AM|^^. No- 
tice that these terms are l/Nc suppressed but appear multiplied by large relative 
m\/m^ factors. The combination AM|j^ has no contributions from countert- 
erms of order 1/iVc. This eliminates a potentially large l/N^ uncertainty. 

Therefore, EM corrections to pion and kaon masses have larger uncertainties than 
the combination AM|j^. 

The above numbers are for the current quark masses defined at the scale A = 
1 GeV. There are contributions from QED running of pseudo-Goldstone boson 
masses both in the short- and long-distance counterparts. The short-distance 
ones are discussed in Section Ol and can be obtained from Eqs. ( |3.6|) . They 
are the terms proportional to K^^^ and K^^^ . As discussed previously, there 
is a numerical ambiguity coming from this contribution. This is parametrized 
by the UV scale A in the log dependence of the short- distance counterpart (|3.5|) . 
The reason for the uncertainty is that the mass definition used here corresponds 
to subtracting a QED counterterm corresponding to the integral from A till oo. 
This uncertainty we estimate by varying the scale A. Within the context of the 
ENJL model (see comments in Section |3]l|) , one expects A to vary between the rho 
meson mass and Aenjl, i-e. between 0.8 GeV and 1.2 GeV, roughly. The stability 
region we find for matching between short- and long-distance contributions is in 
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between 0.6 GeV and 0.9 GeV. We take the variation in the contribution of 
and K^^^ by varying the scale A between 0.8 GeV and 1.2 GeV as the uncertainty 
due to the unknown QED counterterm. We estimate therefore this uncertainty 
in the pseudo-Goldstone boson masses to be lower than 1.2 x 10~^ GeV^ for 
individual pion masses and cancelling for (m^+ — m^o)EM- For individual kaon 
masses it is smaller than 7.4 x 10^^ GeV^. For the combination AM^y^ it is 
smaller than 7.5 x 10^^ GeV^, thus negligible. In all cases this is smaller than 
other uncertainties. 

The prediction we get from our calculation for AM^^ due EM virtual correc- 
tions and in the large Nc limit is 

AM|m = (0.98 ± 0.30) ■ 10"^ GeV^ . (4.2) 

We can now add the leading suppressed l/Nc logarithmic contributions. These 
are the log terms proportional to C in [^]. We include them at a the CHPT 
scale u = Mp. Notice that we are neglecting the l/N^ contributions from the 
counterterms. but they cancel, as mentioned above, for the combination AM^y^. 
Using our determination of C from the chiral limit two-point functions the chiral 
logs are the second number below 

: (1.22 ±0.40) ■ 10"^ GeV^ 
: -(0.04 ±0.01) • 10"^ GeV^ 
: (2.32 ± 0.70) • 10"^ GeV^ 
: -(0.010 ±0.003) ■ 10"^ GeV^ 
: (1.06 ±0.32) • 10"^ GeV^ (4.3) 

The combination AM|^ is calculated directly, without using the other results. 
This is why the error is of the same order as for the individual contributions. 
This is the main result of this work. We confirm a large violation of Dashen's 
theorem. Some phenomenological consequences of it are discussed in Section ^ 

In particular the pion mass difference result should be compared with the 
experimental mass difference — m^o = 1.26 -10"^ GeV^. As expected the 
experimental mass difference value — ?tj,^o is mostly saturated by QED virtual 
contributions, with 30% uncertainty though. The uncertainty here due to not 
included 1/Nc suppressed counterterm contributions, is however larger due to 
the mj^Ks term. An estimate of its contribution will be discussed in the next 
section. 



ml+ 


EM 


(0.95 ±0.27) ■ 


10- 


3 GeV' 


mlo 


EM 


(0.00 -0.04) ■ 
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3 GeV' 


m\+ 


EM 


(1.93 ±0.39) ■ 
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EM 


(-0.01 ±0.00) 
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-3 GeV' 
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(0.98 ±0.08) • 


10- 


3 GeV' 



4.2 Determination of Couplings of the e p Lagrangian 

In this section we give the value of the large couplings in (|2.9| ) and (|2.12| ) 
that can be determined from our calculation. Essentially we have fitted the 
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Figure 6: The long-distance (LD), short-distance (SD) and the sum of the 
Nf. ^ oo results as a function of the matching scale /i for Kiq^u = Mp). 



CHPT large N^. results in Eq. ( |A.1| ) to the output of our calculation. Being off 



mass-shell has allowed to determine one more coupling. See the comment about 



the gauge dependence of these couplings in Section pT3| . We give the central 
values at the points where the stability is best and the errors include typical 
error estimates as well as matching uncertainties. See the comment on 



matching in Section |3.2.3[ As customary we take the CHPT scale v to be the 



rho mass. The QED short-distance contribution in ( |3.6| ) is taken at A = 1 
GeV. The matching scale is always between 0.7 GeV and 0.9 GeV. For some 
combinations of couplings we do not get a very good matching contrary to the 
masses themselves. The variations for /i between 0.7 GeV and 0.9 GeV are still 
within 10% for most. This seems to be caused by the large role played by the QED 
mass renormalization effects. As an example we have plotted the long distance, 
the short distance and the sum for the coefficient Kiq in Fig. |^ and a combination 
with very large cancellations in Fig. ^. The other cases are somewhere in between. 

From the charged combinations in the chiral limit, i.e. the zero mass charged 
pion case, we obtain a good matching and we obtain 

= [(0.00 + 0.30 + 0.67 = 0.97) ± 0.30] ■ lO^^ GeV^ 

at /i = 0.85 GeV, (4.4) 



ml, 



where the first figure is the short-distance QED contribution, the second the rest 
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Figure 7: The long-distance (LD), short-distance (SD) and the sum 
of the Nc ^ oo results as a function of the matching scale /i for 
(-2(2^3 + ^4) + 5(^5 + Ke) - 10(^9 + K^) - ISKn) {ly = M^). The units 
are lO"^. 
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of the short-distance and the third is the long-distance part. The difference with 
the charged pion case is very smalL Using consistently the chiral limit ENJL 
valuei Fq = 89 MeV, we obtain for the e^p*^ CHPT scale independent coefficient 

C = [(0.00 + 1.3 + 2.9 = 4.2) ± 1.5] ■ 10"^ GeV^ 

at fx = 0.85 GeV. (4.5) 

The central value differs by about 30% from earlier determinations, e.g. [|T5[] . 
There it was assumed that the full measured pion mass difference came from the 
e^p*^ term proportional to C. Taking into account our error bars both results are 
nicely compatible, especially since the order e^p^ contribution is about 25% - see 
[|1|, m] and the results after Eq. (U). 

Going off-shell we can get one more combination of couplings of order e^p^, 
namely 

(-2(2^3 + K,) + 5(7^5 + Ke) - 10{K, + K,o) - 18Ku) {v = M,) = 
[(0.85 - 2.53 + 2.64 = 0.96) ± 0.4] ■ 10^^ 

at /i = 0.8 GeV (4.6) 

The three first figures are as in ( [4.4| ). 

From the neutral combinations in the chiral limit we obtain another combi- 
nation of couplings, namely 

(2k^ + K^ + 2(7^5 + Ke) - A{Kg + K^o)) {u = M^) = 
- [(0.68 - 0.30 + 0.13 = 0.51) ± 0.3] ■ 10^^ 

at /i = 0.7 GeV. (4.7) 

Then, by including non-zero current quark masses we can get one more cou- 
pling from the neutral combinations, namely 

{Kg + Kio) {u = Mp) = [(1.7 + 0.4 + 0.6 = 2.7) ± 1.0] ■ lO^^ 

at /i = 0.7 GeV (4.8) 

and combining the charged combination with the chiral limit and the neutral 

Kio{u = Mp) = [(1.7 + 0.5 + 1.8 = 4.0) ± 1.5] • 10^^ 

at /i = 0.7 GeV (4.9) 



^The chiral limit value of is determined in CHPT to be Fq = (86 ± 10) MeV. Since our 
low-energy calculation has been done within the ENJL model we use here the value Fq obtained 
in an ENJL fit to several low energy observables |^ . Notice that this value is also compatible 
with the CHPT value. 
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from the charged kaon mass EM corrections. We could have used the Kiq combi- 
nation from the pion mass EM corrections but the errors due to the subtraction 
of the chiral hmit are larger. We use it as a check. Combinations of these four 
couplings in ( [4.6|) -( ^^) can now be used for other predictions provided the same 
gauge and MS scheme we use are used too. 

From the combinations of couplings above, we can obtain two combinations 
which are free of QED uncertainties, namely, 2K3 + K4 + 2(7^^5 + Kq) and Kg. 
In general, we see a strong dependence on the logarithmically divergent short- 
distance QED contribution. This makes more relevant the danger pointed out in 
Section p.l| of making a naive chiral power counting here. Fortunately, as we have 
seen numerically in the previous section, when combined with the mass factors 
and electric charges, this short- distance QED contribution gives very small final 
contribution to the EM mass corrections. 

The combinations above are determined directly from simple combinations of 
our results. We can combine the numbers in (^]6|) to ( |4.9| ) to obtain as well 

Kg{u = Mp) = -1.3■10"^ (4.10) 
Kio{i^ = Mp) = 4.0■10-^ (4.11) 
2K3 + K^ + 2{K^ + Ke)]{i^ = Mp) = 5.7 ■ 10-^ (4.12) 

~2K3 + Ki + 4:Kii\{u = Mp) = -5.0■10"^ (4.13) 

In view of the results we get for the couplings above, neither short-distance 
or long distance alone dominate any of the couplings. There are in fact large 
cancellations in some cases. So that not much can be said about the couplings 
we don't get from our calculation. 



5 Comparison with Earlier Work 

Historically, the soft pion limit was used in the first attempts for estimating the 
EM contributions to the pions. So that, the authors in |jl| arrived to the following 
expression for the charged pseudo-Goldstone bosons. 
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for the pions and exchanging the Gell-Mann's flavour SU(3) matrices X^^'^ / \/2 
by [Ae + {—)i\-j\/2 for the kaons. Neutral pseudo-Goldstone bosons get zero 
contribution in this limit. The sum rule in ( |5.1| ) has a good high energy behaviour 
in QCD due to the fulfilling of the Weinberg sum rules (WSRs) [^. In fact, in 



the chiral limit, the first and second WSRs guarantee its convergency |I|]. To 
lowest order in CHPT the integrands in (|5.1|) for kaons and pions are the same, 
this is again Dashen's theorem. 

Attempts to go beyond the approximation in (|5] 



11, 13, 20, 21, 23, 24 



are m 

Let us compare our results with the ones obtained in these references. 

The chiral logs at z/ = 1 GeV were used as an estimate of the order of mag- 
nitude of the virtual EM corrections to pseudoscalar Goldstone boson masses in 



[ p^ . This is of course a scale dependent statement, and a conclusion at any par- 
ticular scale is a dangerous one. Only after adding the counterterms the result 
makes sense. The present work is devoted to estimate them. 

In our final result ([4.3|) , we get a relatively large violation of Dashen's theorem, 
though not as large as in and the same as in |^. However, as noticed in 
, the calculation in |^ does not have the correct chiral symmetry behaviour. 
The short-distance contribution is here also assumed to be negligible. We find 
that this could be the case for scales larger than a few GeV. There is also work 
4l[] improving the estimate made in ||20| . The results found there are compatible 



211 some VMD-like estimate in the same line as in 20 is done. The 



with ours. 
In Ref. 

matching scale in the resonant saturation in |21| is identified wrongly with the 
CHPT scale, in our notation here this is setting ^ = v. This is very dangerous, 
since the chiral logs have in both effective theories completely different depen- 
dence. Also no attempt to make any matching of the resonant contribution with 
the short-distance contribution was done. These two points make it very hard 
for us to understand the meaning of their final result. 

Splitting the different contributions to the corrections of Dashen's lowest order 
result, we get 
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0.98 + 0.08 
L26 



0.84 ±0.25 



(5.3) 



where the CHPT scale of the chiral logs is taken at u = Mp. The numbers are 
the leading A^^ correction and the 1/A"c suppressed chiral logarithms. 

Let us now compare this with what we would have gotten from a e^p^ calcu- 



lation using the central values of the Ki counterterms determined in Section ^]2 



at A^. 



oo: 



[0.94(7^:) -0.16(77r) + 0.53(A:i) -0.31(^5^)] ■ 10"^ GeV^ 
1.00 ■ lO^^GeV^ (5.4) 
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The four contributions above are the photon loop contribution to the pion and 
the kaon, the contribution from the Ki counterterms and the L^C contribution, 
respectively. This should be compared with 0.98 ■ 10~^ GeV^ from the full calcu- 
lation. The agreement is very good, remember that our nominally e^p^ couplings 
include also corrections to all orders in CHPT. Notice also the large cancellation 
between the and the L^C contributions, as mentioned before. The latter shows 
the danger of including corrections of the known constants and fully neglect the 
others. The statement is of course somewhat scale dependent. That the chiral 
logs give a sizeable part of the total result at u = Mp is a non-trivial dynamical 
statement which our calculation answers. 

In [|TT| , [13| the order chiral logs together with some order of magnitude 



estimate of the counterterms was used. The large scale dependence of the logs 
only allowed to give a very broad range of results, of course compatible with the 
large Dashen's violation we get. 

The main difference with the calculation in p4[ is the inclusion of the 1/Nc 
suppressed logarithms and the term proportional to L5 (last term in ( |5.4| )). 

Recently there has been some lattice QCD results on the EM contributions to 
hadrons using some quenched unimproved Wilson action. Their final result is 
AM|]yi = 0.64-10"^ GeV^. Unfortunately only an estimate of part of the system- 
atic errors due to finite size effects is given there making difficult the assessment 



on 



of the result. This is particularly relevant after the recent re-analyses ||^ 
the lattice QCD light-quark masses obtained using improved Wilson actions to 
calculate the same observables as in 0. Large 0{a) lattice spacing effects have 
been reported in those works. 

There have been several other calculations in the NJL model done. These 
were all performed at Gy = 0, keeping only the scalar four-quark operators in 
( |3.16|) . One does however expect already at scales around 0.5 GeV (as we 
observe) non-negligible contributions from spin one structures. In addition they 
only treated the pion case, therefore we have not done a full comparison of our 
results with those in 



We can compare the determination of the couplings from the previous section 
with the estimates made in and [^, |T3|. We get that the contribution of 
the counterterms to AM|]y[ is roughly one order of magnitude smaller than the 
estimate made in those references. 

Though the coupling is l/N^ suppressed, its contribution to the charged 
pion mass is potentially large due to a factor m\/m^, as mentioned before. 



We can estimate the short distance estimate to f^g from ( |3.10|) to be around 
Kf^iv = Mp) = -(1.0 ± 1.4) ■ 10"^. An estimate of its total value can be 
done using the result in Table |1] for the pion mass difference and assuming that 
the deviation from the experimental result is just due to the l/N^ counterterms 
proportional to rrij^. Notice that we can do this because we have the complete 
leading in l/N^ contributions. In that way we get Kgi^u = Mp) = —(0.8 ± 2.0) ■ 
10"'^, both compatible with the short- distance estimate and Zweig's rule. In 
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p^l this same coupling was estimated assuming that this deviation is dominated 
completely by the order e^p^ Ki counterterms contributions (both leading and 
next-to-leading in l/Nc) proportional to m|-, so that they get Kgi^i' = Mp) = 
-(4.0 ± 1.7) ■ 10-\ 



6 Ratios of Light- Quark Masses 



A combination of light current quark masses without the Kaplan-Manohar 
ambiguity, can be obtained using the following relation 0, 



m 



+ mlo 




K+ 



2{m 



m 



QCD 



Ami ~ {f^u + 
4 {ml - ml) 



1 



0{ml/k 



(6.1) 



As emphasized in the higher order corrections to this particular ratio are 
very suppressed, so it is very constraining. 

In previous sections, we have computed the virtual EM corrections to the 
pseudo-Goldstone bosons masses, the results are in Eq. ([4.3|) . Subtracting them 
from the experimental masses we get what would be the QCD values for those 
masses. So the result we get for Q is 



Q = 22.0 ±0.6 



(6.2) 



We have calculated the EM corrections at leading order in l/iV^ and to all orders 
in CHPT. The long distance contributions are estimated in the ENJL model and 
the short- distance contributions in the large A^^^ limit of QCD. We also add the 
1/A^c suppressed chiral logarithms to order e^p^. The uncertainty corresponds to 
a 30% uncertainty on our estimate of AM^j^. This should be compared with the 
lowest order result 

Q = 24.2 (6.3) 

using Dashen's theorem. 

The 7] — >■ 7r+7r~7r° decay rate is inversely proportional to the fourth power 
of Q m Recently there has been some activity in improving the prediction 
of the proportionality factor Comparing their best estimate with the 

experimental data, they get 



Q 
Q 



22.4±0.9[||] 



22.7 ±0.8 47 



(6.4) 



in very good agreement with our result. 

To obtain ratios of quark masses themselves requires extra information beyond 
CHPT. A very recent discussion on the consequences of values of Q for the ratios 
of the light current quark masses can be found in EH . 
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7 Polarizabilities of Pseudo-Goldstone Bosons 



In a previous Letter we calculated the cross-section for 77 — >■ Tr^TT*^ to all orders 
in the chiral expansion and leading order in 1/Nc within the ENJL model [Q. 
As part of the work needed for this paper we need to calculate the pseudoscalar- 
pseudoscalar-vector-vector (PPVV) four-point function within the same model 
also for charged pseudoscalars and any flavour channel. We can then obtain 
predictions for 77 7t~^7i~ , K~^K~ , K^K^ to all orders in the chiral expansion 
and leading in the l/N^ expansion. The cross-section for the 77 7r~^7c~ process 
is however dominated by the Born term and higher corrections to it are small [0, 
49[| . The possibility of extracting the pion polarizabilities from the 77 PP^ 



cross-section is studied in [^]. Combining CHPT with dispersive methods they 
arrive to the conclusion that the predicted polarizabilities are in agreement with 
the experimental results. They, however, find low sensitivity in the cross-sections 
to the pion polarizabilities (especially to the neutral ones), so that we cannot 
expect a precise determination from the cross-section 77 PPK The processes 
77 — s> K~^K~, K^K^ occur at a too high center of mass energy for the CHPT 
predictions to be reliable. 

The polarizabilities for the lowest pseudoscalar mesons do however fall in the 
regime where we expect CHPT to work. The chiral calculation for the neutral 
and charged [0] pion polarizabilities have been performed to the two- loop level. 
In those works the order counterterms needed were estimated by using the 
resonance saturation model. Though this model has given good results to order 



p'^ in the strong sector |15|, not much is known about its reliability at 0{p^) (and 
higher). It is therefore important to compare its predictions with other models, 
which like the ENJL model we use, also reproduce the success of the resonance 
saturation predictions for the order couplings in the strong sector. Here we 
will provide the ENJL model estimates for them. 

The polarizabilities for both P"*" and P*^ are defined by expanding the Comp- 
ton amplitude in photon momenta near threshold: 



T = 2 



el-e2*(e - 47rmau;i ^2) - 47rm/3 (gl x el) ■ (g^ x e^*) H ; (7.1) 



The phase convention we use can be obtained from this amplitude definition. 
Here m is the pseudo-Goldstone boson mass, and q= {uj,q), e are the photon 
momentum and polarization vector, respectively. In terms of the relativistic 
amplitudes defined in Eq. (2) of [n] (we use the same notation as there), the 
polarizabilities for both neutral and charged pseudoscalar bosons are given by 

a-^ = —\im(A(s,u = s) + 8m^B(s,u = s)) 

a + (3 = —\imm^B(s,iy = s). (7.2) 

m s^o 



The barred amplitudes in ( [7.2|) are the corresponding amplitudes with the Born 
contributions using pseudoscalar propagators to all orders in CHPT subtracted. 
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We have numerically calculated the full PPVV four-point function, reduced 
the external pseudoscalar legs (see Ref. p5|) and extracted the A(s, z/) and 
B{s, v) form factors. This was done in the same way as was done in [|1^ for the 
neutral case but now we optimized the extraction for s ~ and t k, . To order 

chiral symmetry imposes that -B(s, v) = 0. To the same order A{s, v) also 
vanishes at large Nc for the neutral pseudo-Goldstone bosons (They have 

however a non-zero 0{p'^) contribution from chiral logs |^9|.) For the charged 
cases (tt"^ and K^) we get at large N^., i.e. only from the counterterms, 

A'^''\s,u) = l.l?, GeV-2 (7.3) 

at leading order in the l/iV^ expansion. This compares well with the recent 
determination ^ 

A(^)(s, 1^) = ^ {Lg + Lio) = 1.5 ± 0.2 GeV-^ (7.4) 

TT 

CHPT predicts the following counterterm structure for the form factors A{s, u) 
and B{s, v) 

A{s, u)count. = A^-^\s, u) + . + ■ ■ ■ , 

b 

Bis, v)count. = 7-^774 + " " " (7-5) 

where Fi is for pions and for kaons. In the ENJL model Fq = 89 MeV, 
= 90 MeV and F^ = 96 MeV at large N^. The second term in A{s, v) 
and first in B{s,v) are the order contributions. We have fitted the A{s,v) 
and B[s, v) ENJL form factors to a polynomial in s up to order s^. These give 
always reasonable good fits. Thus, the z/^ dependence in the energy region where 
polarizabilities are defined is small. The results of those fits to the form factors 
at A{s, v) and -B(s, v) leading order in l/N^ in terms of ai, 02 and bi are given 
in Table |^ (here we consistently used the ENJL model values for Fq, F^ and F^ 
given above). We have done it for the chiral limit pseudo-Goldstone boson, the 
pion and the kaon with their appropriate masses, both for the charged and for 
the neutral case. The results we get from the fit for the coefficients ai, 02 and 
bi in Table ^ include higher than order p^ corrections which are not of the type 
s^, i.e. mainly mass corrections. The coefficients in Table have a typical error 
estimate for the l/N^ expansion of (20 ~ 30) %. In fact the difference between 
the coefficients 02 and bi for the vr"*" and the are higher order corrections. The 
coefficients ai of the charged pion and kaon have larger uncertainty than the rest 
since we get them from subtracting the dominant order p^ contribution from the 
A{s, v) form factor. From the coefficients of the chiral limit pseudo-Goldstone 
boson in Table ^ and the ai coefficients for the charged pion and kaon, one can 
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ai 


a2 


bi 






6.7 


0.38 


71+ 


-8.7 


5.9 


0.38 




-5.6 


15.8 


0.77 


x' 




14.0 


1.66 




-23.3 


14.9 


1.69 




-13.2 


16.9 


1.10 



Table 2: The dimensionless ai, 02 and bi coefficients for the chiral hmit 
pseudo-Goldstone boson x, the pion and kaon, both for the charged and neu- 
tral ones. 





and Higher Count. 


/ X-logs 


Total 




0.62 


0.17 


0.79 ± 0.25 


7r+ 


0.14 


0.31 


0.45 ± 0.15 




1.16 






K+ 


0.81 







Table 3: The combination of polarizabilities a + (3 in units of 10 ^ fm^. The 
chiral logs for the kaons (Column 3) are not known. 

obtain the ENJL predictions for the six terms of the order chiral Lagrangian 
p5[] contributing to 77 PP"'' at large Nc- 

To get the complete prediction for the polarizabilities, one has to add to the 
counterterm contribution we calculate within the the ENJL model to all orders 
in the CHPT expansion, the contributions from chiral loop diagrams (order 1/Nc 
in the large Nc counting). These are for the pions known to two loops P, |^. 
Our final result for the pion and kaons polarizabilities are in Tables |^ and ^. For 
the SU(2) Ij counterterms (see for their definition) entering the pion chiral log 
expressions we have consistently used the ENJL predicted ones, i.e. 

Ii = -1.05, I2 = 5.8, I3 = 2.5, I4 = 4.3, I5 = 14.6, Ig = 16.9 . (7.6) 



The final result in Tables |^ and ^ contain the counterterm contributions to 
all orders in CHPT and large Nc. Notice that the experimental result from (|7.4| ) 
is {a — l3)o{p*') = (6.0 ± 0.8) ■ 10~^ fm^. We can compare the counterterm entries 
in Tables § and § with results obtained using resonance saturation and other 



calculations in the ENJL modelfSlJ. For the comparison between the resonance 



saturation estimates of for the neutral pion case see Table 1 and text in [|4 
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p'^ Count. 


X-logs 


and Higher Count. 


/ X-logs 


Total 







-1.08 


-0.47 


-0.31 


-1.9 ± 0.6 


7r+ 


4.53 





-0.25 


-1.13 


3.2 ± 1.1 


i^o 







-0.58 






K+ 


1.27 




0.07 







Table 4: The combination of polarizabilities a — (3 in units of 10 ^ iw?. The 
chiral logs for the kaons (Columns 3 and 5) are not known. 

In 1^ the following results where obtained from resonance saturation for the 
coefficients in Table @ in the case of the vr"*". 

ai = -3.7 ±1.65; as = 0.75 ± 0.65 ; 6i = 0.45 ± 0.15 . (7.7) 

Only hi seems to be in agreement, notice though that the estimates in don't 
include the contributions from scalar and tensor resonances. The main part of 
that work is however the two-loop calculation of the charged pion polarizabilities. 
We disagree, as discussed in , with the way the order coefficients were 



obtained in [^]. We have not compared with the calculations in |5^, ^ because 
they work in the ENJL model with Gy = 0, so the important (even dominant) 
effects coming from vector and axial-vector mesons exchanges are not included. 
For comparison with predictions of other models for the pion polarizabilities see 

For a recent review of the experimental situation and data on pion polariz- 
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abilities see R 54 1 and references therein. 



8 Summary and Conclusions 

The main aim of this work has been the calculation of the virtual EM corrections 
to the masses of the pseudo-Goldstone bosons. This was motivated by some re- 
cent calculations pO|, |2^, Q where large corrections to Dashen's theorem were 
obtained. This was also supported by recent improved calculations of the decay 
rate for r] tt+tt^vt^ p6| , |47[| which included estimates of higher order corrections 
using dispersive techniques. At the same time there appeared some works claim- 



ing that small violations of Dashen's theorem were not excluded [11, 13, 21]. We 



have calculated in the large Nc limit and using a technique similar to the one in 



1^ but for Green functions off-shell, the virtual EM mass corrections to pions 
and kaons, setting rriu = rrid- Our result is 

= (1.06 ± 0.32) • 10"^ GeV^ (8.1) 
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where we have included the known 1/Nc suppressed chiral logs at order e^p^ 
(they are only 0.08 of that number). The error includes an estimate of the 1/iVc 
corrections we are missing among other uncertainties discussed along the text. 
Notice that the 1/Nc contribution coming from the order e^p^ counterterms cancel 
in the combination AM|j^. 

Our general conclusion is that a large violation of Dashen's theorem is quite 
well established. At the CHPT scale v = Mp this is dominated by the photon 
loop contribution (both logs and constant pieces). In fact, our calculation is to 
all orders in CHPT (in the long-distance part) at large iVc, so it includes that 
contribution to all orders. The dominance of the photon loop contribution at 
V = Mp is due to a large accidental cancellation between counterterms of order 
e^p^ of both types: proportional to LiC and Ki. The inclusion of just part of 
them is very dangerous at any scale. Variation of the CHPT scale in the logs 
should be accompanied with the running of the counterterms which could become 
eventually important. 

A small remark here. We find some 80 % correction to the EM contribution 
to kaon masses (mainly from the next-to- leading order). This is very similar 
and consistent with what we found in another kaon self-energy quantity, the 



so called i?i^-parameter [^. We observe then that two-point function kaon 
self-energies from gauge-bosons exchange have very large higher order CHPT 
corrections. Notice that in our approach we are able to make a calculation to all 
orders in CHPT at large Nc- 

We have obtained also the ratio of light quark masses defined in ( |6.1| ), 

Q = 22.0 ±0.6, (8.2) 



in good agreement with the one found in Our result supports the very 

recent scheme of ratios of light-quark masses presented by Leutwyler in . 
We have also estimated some couplings of the order e^p^ effective Lagrangian 



described by Urech |TT| . These were all of the expected order of magnitude. For 



values of the combinations obtained here, see Eqs. ( [4. 61) to ( [4.10|) . 



We have discussed the ambiguity of the electromagnetic gauge choice in the 
definition of these couplings and pointed out how in our approach this is cir- 
cumvented. We also discussed how to include the short distance renormalization 
needed due to photon loops. 

Given the large cancellation observed here between the Ki and the Li C coun- 
terterm contributions to electromagnetic mass differences, one should also be 
careful with estimates of the electromagnetic corrections to other quantities when 
only one of these is taken into account. 

As a by-product we have also predicted the order terms which enter in the 
description of 77 PP'^ (P = tt"*", 7r°, 7^°, i^"'") decays and the counterterms of 
electric and magnetic polarizabilities to all orders in CHPT for pions and kaons. 
These predictions were given in Section |^. 
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EM Corrections to Pseudoscalar Two-Point 
Functions to 0{e^p^^ 



In this appendix we give the large expressions for the electromagnetic contri- 
bution to the pseudoscalar two-point functions in (|2.1|) after reducing to order 



p in CHPT. These are finite quantities, the coupling constants Ki and Ki are 
the renormalized finite parts in the MS scheme of |^ at some scale u (this scale 



has nothing to do with the scale /i introduced in (3^) or u in Section |3.1|) . 

We only give here the large Nc expressions off mass-shell, the complete on- 
shell expressions to order p"^ can be found in |ll|] (with the translation of 
couplings in ( |2.15D ). These expressions are the needed ones to obtain the large 
Nc predictions for the Ki and Ki counterterms in ( p.l2| ) from our calculation. We 
also will use the qq basis with q = u,d,s quarks which is more directly related to 
our large Nc calculation. We use the Feynman gauge for the photon propagator. 
The reduced pseudoscalar two-point functions give the following EM corrections 
to the pseudo-Goldstone boson masses 



16e2 
~9~ 



q^ ((2i^3 + 



+ 2{K, + K,)-4{K, + K 



10, 



2ml{K, + K 
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(A.l) 



B Large Nc Short-Distance Contributions 



In this appendix we give the short- distance part of the EM contributions to the 
pseudoscalar two-point functions in in the large N^. hmit. These are the 

contributions of the higher than /i photon modes to the EM corrections to these 
pseudoscalar two-point functions after reducing. (We only give the independent 
ones in the large limit.) 
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